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WIDE AREA COVERAGE RADAR IMAGING SATELLITE FOR EARTH APPLICATIONS 
by Grady H. Stevens and James R. Ramler 

Lewis Research Center 
Cleveland, Ohio 

ABSTRACT 

A preliminary study was made of a radar imaging satellite for earth 
applications. A side-looking synthetic-aperture radar was considered and 
the feasibility of obtaining a wide area coverage to reduce the time re- 
quired to image a given area was investigated. Two basic approaches were 
examined; low altitude sun-synchronous orbits using a mul t i beam/mul t i - 
frequency radar system and equatorial orbits up to near-synchronous altitude 
using a single beam system. Surveillance and mapping of ice on the Great 
Lakes was used as a typical application to focus the study effort. 


ii 



INTRODUCTION 


. The National Aeronautics and Space Administration's continued interest 
in and support of space applications has resulted in recent years in sig- 
nificant advancements in developing technology for earth observations 
including resource surveys, meteorology, geology, oceanography, and environ- 
mental monitoring, all on a global scale. Weather imaging satellites, such 
as TIROS and NIMBUS, have been extremely useful to the advancement of the 
meteorological sciences. The first Earth Resources Technology Satellite 
(ERTS-l) launched on July 23, 1972 is providing impressive data on the 
earth's environment and its resources. Such programs are demonstrating the 
technology required for eventually establishing an operational system of 
earth observation satellites. 

To date, observation satellites have operated in the UV, visible or IR 
spectrums. It is logical to consider complementing these systems with 
observations in the RF spectrum, in particular using radar systems. Some 

a 

activities have already started in this area. For example, a microwave 
experiment operating at about 14 GHz was recently flown on Skylab to measure 
radar scattering cross section and passive microwave emissivity of land and 
sea and to determine signal correlation properties. Satellites in early 
mission definition phases include SEASAT to be launched about 1978 to 
demonstrate and test remote sensing instrumentation for use on oceanographic 
satellites; Earth Observatory Satellite (EOS) to be launched about 1978 as 
a follow-on to ERTS; and SEOS, an experimental synchronous earth observation 
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satellite to be launched about 1S80 with active and passive microwave 
scanners for observing phenomena 9 such as d situations, sea state and 
anomalous snow packs. The LLS. Navy is studying the feasibility of using 
radar imaging satellites to provide sea state and fleet status information. 

While radar has been used for some time for a variety of detection and 
tracking applications,, it has only been recently that its application to 
fine-resolution aerial mapping has been fully developed. The well-known 
capability of airborne radar to generate imagery at night and through cloud 
cover enhances the att ract \ veness of this technique. 

Side-looking radar in which the earth is viewed at oblique angles is 
particularly suited for discerning surface features. It has been developed 
to the point where it is now a proven instrument for all-weather reconnaissance 
and ground mapping. 

Concurrent with the development of side-looking airborne radar has been 
a growing interest in its potential for a wide variety of earth applications. 
These include surveillance and mapping of ice in shipping areas, such as the 
Great Lakes and Alaskan ports, scientific observation of the Arctic ice fields, 
ocean state surve i 1 1 ance, pollution detection and disaster assessment. 

An investigation of the Great Lakes ice cover by the Radar and Optics 
Laboratory of the University of Michigan (ref. 1) using synthetic aperture 
radar suggested that a suitably designed radar system may have the capability 
to discriminate among a large variety of ice types. Such a system could 
potentially impact shipping and ice clearing operations in a major way and 
enhance the economics of the entire region. 
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The Lewis Research Center is currently conducting a demonstration 
program to develop a rapid all-weather ice information system for determining 
ice type and coverage on the Great Lakes for shipping and navigation pur- 
poses. This study is part of a multiagency demonstration program in response 
to Public Law 91-611 passed by Congress in 1970 which authorized a 3-year 
demonstration program for determining the feasibility of extending the Great 
Lakes shipping season. 

During the 1972-73 winter season an AN/AP$-9^C Side-Looking Airborne 
Radar (SLAR) installed in a Lewis Research Center 0V-1B aircraft was flown 
over Lake Erie and the adjacent Great Lakes to obtain radar imagery of the 
ice cover (ref. 2). Results to date have been excellent and it is planned 
to continue this project through the I973“7^ winter season. It is anticipated 
that visual and automatic methods will be developed for interpreting this 
imagery for ice type, distribution, thickness, and motion. As an aid for 
developing these interpretive methods, extensive “ground truth" and low 
altitude visual and thermal imagery are acquired simultaneously with the 
radar imagery. The eventual use of a radar imaging satellite to collect such 
information is a clear possibility. 

Any serious attempt to extend the general shipping season beyond the 
present duration of about 8 months must begin with a modern updated system 
which provides, as often as needed, a map of the ice coverage and ice type 
distribution everywhere on the lakes. At present, only the most primitive 
means are used to provide this information and, unfortunately, no agency 
exists which recognizes as part of its mission, the upgrading of existing 
methods. Such an upgrading requires the type of technology already developed 
and being developed by NASA. 
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This report covers the results of a preliminary feasibility study of a 
radar imaging satellite for earth applications. The purpose of this study 
was to scope the problems associated with radar imaging from orbital 
altitudes, to grossly define the requirements of a radar imaging satellite 
and to assess the technical feasibility of such a system. 

A s ide-1 ook i ng, synthetic aperture radar was considered and the 
feasibility of obtaining a wide area coverage (large swath width) to reduce 
the time required to image a given area was investigated. Two basic 
approaches were examined; low altitude sun-synchronous orbits using a multi- 
beam radar system and equatorial orbits up to near-synchronous altitude to 
reduce the beams required. While such a satellite would likely have a multi- 
application role, surveillance and mapping of ice on the Great Lakes was 
used to focus the study effort. 

Future studies will be needed to identify other potential users of the 
imagery and to determine what applications are practical. The costs and 
benefits associated with these applications will need to be defined along 
with the capabilities of alternate approaches such as aircraft, in situ 
observations, etc. to completely assess the merits of a radar imaging satellite. 
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SUMMARY 

A preliminary study was made of a radar imaging satellite for earth 
appl i cat ions . A s i de- 1 ook i ng , synthetic aperture radar was considered 
and the feasibility of obtaining a wide area coverage (large swath width) 
to reduce the time required to image a given area was investigated. Short 
response time would be necessary for a number of applications, such as 
ice surveillance in shipping areas, ocean state surveillance and pollution 
mon i tor i ng . 

This study was conducted in support of the Lewis Research Center's 
participation in a multiagency demonstration program to determine the 
feasibility of extending the Great Lakes shipping season. Lewis is currently 
flying an AN/APS-94C Side-Looking Airborne Radar (SLAR) over Lake Erie and 
the adjacent Great Lakes to develop a rapid all-weather ice information 
system for determining ice type and coverage on the Great Lakes for shipping 
and navigation purposes. Consequently, surveillance and mapping of ice 
on the Great Lakes was used as a typical application to focus the satellite 
study effort. For ice imaging the conditions being surveyed can be quite 
changeable due to current, winds, temperature, etc. and it would be 
necessary to obtain a complete image in as short a time as possible (daily 
if possible) to effectively aid navigation and ice control operations. 
However, it is obvious that this single application would not justify a 
dedicated satellite and that it would likely perform a multiapplication 
role . 

The purpose of this study was to scope the problems associated with 
wideswath radar imaging from orbital altitudes, to grossly define the 
requirements of such a satellite and to assess its technical feasibility. 
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Two basic approaches were examined; low altitude sun-synchronous orbits 
using a mul t i beam/mul t i frequency radar system and equatorial orbits up to 
near-synchronous altitude to reduce the number of beams required. Equatorial 
orbits do not appear to be a practical solution. Two equatorial orbits were 
considered in some detail; one at about 6400 km altitude and another at 
about 14,000 km. 

The 6400 km orbit "laps" the earth five times daily thus allowing a 
daily image to be constructed from a mosaic of five consecutive passes. 

Radar beam swath widths of about 250 km are achievable at the latitude of 
the Great Lakes with such an equatorial satellite so that an area the size 
of the Great Lakes (about 900 km north to south) could be mapped in five or 
even four passes. For the single application of mapping the Great Lakes 
(or any area at the same latitude) a single beam/single frequency synthetic 
aperture radar system would be required with an average power requirement of 
about 400 watts . 

However, an equatorial orbit at 6400 km altitude would place the satellite 
in the midst of the most intense trapped proton radiation belt and a power 
source other than solar cells (e.g., an RTG) would be required. Other 
electronic components would need radiation protection, thus incurring a 
weight penalty. This system is judged to be relatively complex and costly 
and does not appear to be an attractive solution. However, this approach 
would need more detailed study than can be given in this report to adequately 
assess its feasibility. 

An alternate approach was studied in which the orbit altitude would be 
high enough so that solar cells could be adequately shielded. An altitude 
of about 14,000 km was chosen. At this altitude, daily imaging with a single 
beam/single frequency radar is no longer possible without resorting to 
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antennas in excess of 30 meters. The required average power for imaging 
at the latitude of the Great Lakes would be increased to about 2.3 kw. 

The satellite would "lap 1 ' the earth twice daily. At the Great Lakes, swath 
widths of 280-420 km appear feasible depending on antenna size. It would 
take from 3-4 days to map an area the size of the Great Lakes with the 
mosa i c technique. Battery backup power does not appear practical for 
eclipse periods due to the high daily energy demand. 

Equatorial orbits do not allow global coverage thus limiting other 
potential applications. Coverage of a range of latitudes would require a 
mu 1 1 i beam/mu 1 1 i frequency synthetic aperture radar and a slewed antenna. In 
addition, coverage at the lower latitudes results in high radar beam 
incidence angles. In order to maintain a constant range resolution as the 
latitude is decreased the power would need to , be increased markedly up to 
the mu 1 1 1 k i 1 owatt range. Thus, the lowest latitude covered would determine 
the power requirement. 

Sun-synchronous orbits around 1000 km altitude appear to be the best 
approach to obtaining wide swath coverage. However, a mul t i beam/mul t i frequency 
synthetic aperture radar system will need to be developed. Daily coverage 
of an area the size of the Great Lakes does not appear practical since the 
large number of beams required (approximately 12) is not consistent with 
the required antenna range dimension of about 3 meters. Coverage at 2-day 
or 3 _ day intervals appears feasible. This would require six beamsand four 
beams respectively. The average power required would be about 600 and 400 
watts respectively. 
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The 3-day and 2-day coverage cases are summarized below for an X-band 
radar with a resolution of 100 meters x 100 meters and a receiver effective 
noise temperature of 500 D K. 

Sun-Synchronous Orbits 


Coverage Period 

3-day 

2-day 

Altitude, km 


101 1 .5441 

1072.2403 

Period, hours 


1 .75610 

1.77778 

1 ncl i nat i on, deg 


99.5310 

99.8105 

Average power, watts 


400 

600 

Orbi ts per 24 hours 


13-2/3 

13-1/2 

Orbit track separation 

in longitude, deg 

8.7805 

13.3333 

Orbit track separation 

at equator, km 

952.0 

1443.5 

Orbit track separation 
km 

at 45 deg 1 at i tude, 

655.5 

995.2 

Beams required 


4 

6 

Total achievable swath 

width, km 

436 

654 

Antenna size, length x 

width, meters 

10x3 

10x3 


Sun-synchronous orbits can obtain global coverage, thus allowing 
potentially more applications than an equatorial orbit. Antenna slewing 
would not be required since all targets eventually come into view. 

Future studies are needed to define this system in more depth. Par- 
ticular attention should be given to assessing the feasibility of a multibeam/ 
mul t i frequency synthetic aperture radar and the possibility of using a large 
number of beams such as 12 to allow daily coverage of an area the size of 
the Great Lakes. In addition, studies are needed to identify other potential 
users and to determine what applications are practical. The costs and 
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benefits associated with these applications need to be defined along with 
the capabilities of alternate approaches such as aircraft, in situ observa- 
tion, etc. to adequately assess the merits of a radar imaging satellite. 
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ORBIT CONSIDERATIONS 
A. General r 

Radar Imaging from a satellite Is basically a different and in many 
ways technically more demanding application than using aircraft mounted 
radar. Because of the higher altitude, considerably higher RF power levels 
are required. In addition, some satellite orbits result in a much higher 
ground velocity than an aircraft which can severely limit the maximum ground 
area covered on any one pass for a given set of parameters such as image 
resolution, RF power, etc. 

Initial rough estimates of the radiated RF power requ i red indicated 
that it could easily be on the order of several kilowatts so that it is 
important that a reasonably optimum system be devised to limit the power 
required to a practically achievable value. In general, the lower the 
altitude the lower will be the power required (all other things being equal). 
This suggested that one class of orbits to be considered should be low earth 
orbits and in particular, sun-synchronous orbits. Sun-synchronous orbits 
allow the satellite to remain in more or less continuous sunlight so that 
the relatively high power levels required can be obtained from solar cells 
on a nearly continuous basis. One major problem that low sun-synchronous 
orbits have, however, is that the velocity of the satellite (or radar beam) 
relative to the earth's surface is quite high (approximately orbital velocity) 
This severely limits the ability of the radar system to obtain a good image 
resolution unless the radar beam is restricted to a fairly small swath as 
it sweeps across the earth's surface. This is discussed in detail in the 
section dealing with Radar Considerations. This swath could be extended by 
using a multi-beam system with an associated increase in the RF power required 
However, such a system may be complex and costly to implement. 



One possible means for avoiding these problems is to reduce the 
satellite's velocity relative to the earth's surface thereby improving 
the ability of the radar system to cover a much wider swath while retaining 
the desired image resolution. This would suggest relatively high equatorial 
orbits up to near-geosynchronous altitude where the earth-relative velocity 
could be tailored to some desired value, (A geosynchronous equatorial orbit 
has no earth-relative velocity and could not be used since the radar depends 
on a Doppler effect for resolution.) As the orbit altitude is increased 
a wide swath width can be obtained with fewer beams or even a single beam 
system thus reducing the power required as the number of beams is reduced. 
However, at the same time the power will increase as the altitude is 
increased. It then becomes a matter of trading off these effects to find 
an optimum equatorial orbit. 

The following sections will discuss some of the factors to be con- 
sidered for the two basic types of orbits that seem feasible: (1) low sun- 

synchronous and (2) equatorial. 

B. Sun-Synchronous Orbits 

A sun-synchronous orbit is defined as an orbit which has a nodal pre- 
cession rate equal in magnitude and direction to the earth's mean rate of 
revolution about the sun, 1.991 X 1 0 ~ 7 rad/sec (O .9856 deg/day) „ The earth- 
sun-satellite geometry is shown in Figure 1. This figure depicts a geo- 
centric coordinate system in which the X-axis is directed toward the vernal 
equinox, the Z-axis is directed toward the north celestial pole, and the 
Y-axis lies in the equatorial plane forming a right-handed coordinate system 
with the X and Z axes. The orientation of the orbit plane is described by: 
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a The node ( ) in the equatorial plane measured eastward from the 

vernal equinox to the ascending node of the orbit plane, 
b. The inclination (1) measured in a counterclockwise direction at 
the ascending node from the equatorial plane to the orbit plane. 

The ecliptic plane is the plane of the motion of the sun as seen from the 

center of the earth. The obliquity of the ecliptic (€) is the angle 

between the ecliptic and equatorial planes. The position of the sun in the 
ecliptic plane is described by the celestial longitude of the sun (X) which 
is measured eastward from the vernal equinox to the sun. The position of 
the satellite in the orbit plane is described by the argument of latitude (u) 
which is measured from the ascending node to the satellite in the direction 

of motion. 

The desired nodal precession of the satellite orbit can be accomplished 
by utilizing the effects of earth's oblateness. An analytical expression 
for the secular rate of change of nodes of an orbit can be obtained by con- 
sidering only the terms through the second harmonic in the potential function 

of the oblate earth. The following expression is obtained: 

& -/p" JR e 2 ( R e + h )' 7/Z cos l for 0 Si S 2 * (1) 

where 

3 

u = gravitational constant of earth, 398601 .2 knri_^ 

J = coefficient of second harmonic in gravitational potential function, 

1 .62405 X 10“ 3 

R e = earth equatorial radius, 6378.160 km 
h = altitude of satellite above equatorial radius 


1 = orbit inclination 
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For inclinations greater than jt rad (retrograde orbits), A , is positive, 

2 3 

which results in a precession (eastward movement) of the node. This is the 

desired direction to keep the orbit plane in a proper position relative to 

the sun. Inclinations less than jx (posigrade orbits) will result in a re- 

2 

gression (westward movement) of the orbit plane. Perturbations due to the 
oblateness of the earth on other orbital elements such as inclination, 
eccentricity, and semimajor axis are small and periodic with no resultant 
secular change. 

A plot of equation (1) is shown in Figure 2 for the desired precession 
rate of 1.991 X 10 ^ rad/sec. Sun-synchronous orbits about the earth do 
not exist above 597^.7 km (3226.1 n mi). It is obvious that for relatively 
low earth orbits (altitudes less than 2000 km say) the orbits will be nearly 
pol ar . 

Consider the geometry (Figure 3) of a satellite in a nearly polar 
retrograde circular orbit about the earth. THe orbital angular velocity 
vector is defined as the orbit perpendicular. The longitudinal position or 
right ascension of the orbit perpendicular ( ) is measured positively 

counterclockwise from the vernal equinox to the projection of the orbit 
perpendicular in the equatorial plane. Note that + 3^/2. The 

sun (sunline) may in general be on either side of the orbit plane. For this 
analysis it was taken to be on the same side as shown in Figure 3. The 
opposite side analysis is similar. The orientation of the sunline is described 
by angles obtainable from an ephemeris: 

a. The right ascension of the sun (£<, ) measured positively counter- 
clockwise from the vernal equinox to the projection of the sunline 
in the equatorial plane. 
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b. The declination of the sun (<n measured positively above the 
equator and negative below. 

The longitudinal relation between the sunline and the orbit perpendicular 
can be taken as ^ ~ Ot - • A so-called noon-midnight sun-synchronous 

orbit is one in which ^ is nominally n/% or 3^/2 • Some sun-synchronous 
satellites such as the ITOS meteorological satellite, use an orbit in which 
the satellite always crosses the equator at 3p.m. northbound and 3a.m. south 
bound local time ( ^ = approximate! y Jt/4) . This type of orbit provides 
oblique illumination angles suitable for photographic purposes. 

In the case of a radar imaging satel I i te, i 1 1 um J nat ion of the earth is 
not required and the orbit may be chosen to provide maximum illumination of 
the satellite's solar cells. Twilight orbits where V* is nominally zero 
are ideally suited for this purpose since the satellite would be illuminated 
over all or almost all of the orbit. Any desired V“* may be established by 
simply launching at the proper time of day. 

1 . Orb it Eel i pse 

By assuming only circular orbits and a spherical earth, the geometry of 
the orbit with respect to illumination from the sun is simplified. The 
existence of a penumbra region of shadow can be neglected with negligible 
error thus assuming that the rays of the sun are parallel, causing a cylinde 
shaped umbra region behind the earth. Shown in Figure 4 is an edgewise view 
of the orbit plane and its geometrical relation to the position of the sun. 
The angle between the orbit perpendicular and the sunline is rj. It can be 
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seen from Figure 4 that as long as rj is less than the critical value, 

the entire orbit will be in sunlight. The situation depicted in Figure 4 

is a limiting case in which the orbit is just on the edge of the cylindrical 

shadow and r\ equals rj . From Figure 4, t \ is defined by 

c c 


cos 7] = for 0 < rj < £ 

c + h c 2 


( 2 ) 


This means that for an orbit to be in continuous sunlight, the orbit per- 
pendicular must remain inside a cone having the moving sunline as its axis 
(ref. 3). The half-angle of this cone is the corresponding to the 
instantaneous orbit altitude. An expression for tj can be derived from 
Figure 3. For the sunline being on the same side of the orbit plane as the 
orbit perpendicular, it can be shown that 

cos T) = cos V* sin 1 cos & + cos 1 sin £ (3) 

The angles rj and r\ can then be evaluated from equations (2) and (3), 
respectively, and compared to determine whether or not the entire satellite 
orb it is in sunl i ght . 

If the orbit is eclipsed, the fraction of the orbit in shadow can be 


simply found from the equation_ 

f = £ + s i n 
2 


-1 



for T] > TJ 


po 


7t 

However, in order to determine battery requirements, charge/discharge rates 
and duty cycles it is also necessary to know when or at what point in the 
orbit the satellite will enter and leave the eclipsed segment. From reference 
4 it can be shown that 
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cos (n/2 + q ) = sin u b in X (sin € fein I + cos € cos I cos ■a a ) 

- cos X (cos 1 sin ■a a] + cos u £sin X (cos £ sinilj a ) + 

cos X (cos ^ a )j (5) 

where sin X = sin <£ 
sin € 

cos X = cos <S - sin A sin c( cos 6 

cos o< 

Two solutions for u may be obtained from these equations corresponding to 
the satellite orbital positions at entrance to and exit from the shadow. 

Continuously sunlit orbits only exist over an altitude range from 
about 1500 km (800 n mi) to 3300 km (1800 n mi). For other altitudes the 
eclipse time varies depending on the time of year. This is depicted in 
figure 5 for the range of altitudes from 185 km- to 1 296 km (100 n mi to 700 
n mi). Some shadowing occurs when the satellite is over the northern hemi- 
sphere for altitudes below about 280 km (150 n mi). 

It should be noted at this point that if the sun were assumed to be 
on the opposite side of the orbit plane from that depicted in figure 3; the 
shadowing characteristics would be opposite to that shown in figure 5* That 
is, the major eclipsing would occur in the northern hemisphere during the 
winter instead of in the southern hemisphere during the summer. For this 
reason, the "oppos i te-s ide" case is not a desirable geometry for ice 
imaging in the northern hemisphere. 

2, Atmospheric Drag 

a. Orbit Decay : Atmospheric drag causes a satellite to decay or lose 

altitude. For all practical purposes, the orbit remains circular and the 
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inclination remains constant. Atmospheric drag decreases with increasing 
altitude until it is virtually zero at a thousand kilometers or so. How- 
ever, for relatively low orbits (several hundred kilometers), which are 
being considered here, atmospheric drag can be a major perturbation. 

The most accurate estimate of the effects of orbit perturbations on 
the motion of a satellite is obtained by integrating the differential 
equations of motion for an orbiting spacecraft while considering all 
perturbing forces such as atmospheric drag, grav i tat ional perturbation of 
other bodies, solar pressure etc., as precisely as possible. However, this 
requires extensive input data and computer time. A number of techniques 
have been developed to yield approximate results making certain simplifying 
assumptions. The accuracy of these other techniques when considering 
relatively low orbits depends primarily on the assumptions made with regard 
to upper atmosphere density and its variation with solar activity. For 
the purposes of this study the effects of solar activity have been averaged. 

Atmospheric drag effects have been estimated assuming an atmosphere 
model based on the 1962 U.S. Standard Atmosphere (ref. 5). This model 
depicts a typical mid-latitude, year-round condition averaged for the range 
of solar activity that occurs between sunspot minimum and maximum. The 
altitude range of interest lies above 200 km. The 1962 U.S. Standard 
Atmosphere model in the altitude range from 90 to 700 km is termed 
speculative. Although the model does not extend beyond 700 km this is 
sufficient for our purposes. Atmospheric density is presented in figure 6 
for an altitude range from 200 to 700 km. 

Following the analysis of Billik (ref. 6) and that of reference 4, 
an exponential atmosphere will be assumed. 
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e = 


fe. 


C.. 

- £ h 


( 6 ) 


where f = 1.58423 X l(f 9 , kg/m 3 
Z 3 = I .33812 X 10' 5 , m' 1 
^ = density, kg/ni 
h = altitude, m 

This assumed fit is shown in figure 6. Note that this is an approximation 
to the actual density but is reasonably close for the altitude range being 
cons i dered . 

For h < < R g , Billik has derived an expression involving h, t and 


known constants 
86400 




(7) 


where W/C^A is the ballistic coefficient involving: W, the satellite weight 

2 

(kgf ) , C^, the satellite drag coefficient (dimensionless), and A (m ), the 
effective satellite area normal to the direction of motion and consistent 
with the definition of Cp. The factor 86400 allows t to be input in days. 


Equation (7) may be solved for h. 


h = 1 In 

-86400 

B / t \ + 

f 


L w/( w J 


( 8 ) 


where B = n \fFX = constant 

A plot of equation (8) is presented in figure 7 for various initial 
altitudes. As an example, for a satellite having a weight of 1000 kgf, an 
area of 1 0m and a drag coefficient of 2.0 (ballistic coefficient equals 50) 
it would take about 320 days to decay from an initial altitude of 500 km 
down to 400 km. 
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To determine the amount of de-synchronization as the satellite decays 
it is useful to consider the longitudinal relation between the sunline and 
the orbit perpendicular ( ^ ) . An equation for ^ in terms of t and known 

t 


constants is derived in Appendix A. A plot of 'K q - f' versus 

“wTcTa 


w/c d a 


for 

is shown in figure 8/various initial sun-synchronous orbits (ho = 300 to 

800 km) o Intermediate altitudes are noted on the curves down to 200 km. 

As an example, consider the effects of decay from ho - 600 km to h « 400 km. 

From figure 8, ^ o 
W/C D‘ 

ballistic coefficient of 50 is assumed, Y* q - T = 1.175 rad (63°3 deg) 
and t = 1550 days (4.25 years). 

As an orbit decays, the eclipse fraction will tend to increase. This 

effect is illustrated in figures 9(a)-9(e) for initial altitudes of 400-800 

2 

km and an assumed ballistic coefficient of 50 kgf/m . The time that the 
satellite would spend in shadow on each orbit is shown as it varies through- 
out the year. The mission is shown starting on January 1 with ^ s 0. 
ff the orbit altitude is maintained (by thrusting to cancel drag) the eclipse 
pattern will simply repeat itself year after year without getting any worse 
(assuming sun-synchronization is also maintained). If no drag compensation 
is provided, the orbit will decay with a consequent increase in the time 
spent in shadow each year. 

2 

For a ballistic coefficient of 50 kgf/m the orbit will completely decay 


^ = 0.0235 rad and t - 31 days 

kgf/m 2 W/C D A kgf/m Z 


If a 


and reenter the earth's atmosphere in less than 5 years for initial altitudes 
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up to slightly beyond 600 km (see figure 7). However, the situation 
rapidly improves beyond this point so that at h Q 85 700 km the altitude would 
decay to only about 675 km after 5 years and the effect on the eclipse 
history would not be significant (figure 9(d)). 

b. Drag Compensation : Based on the considerations just discussed it 

is clear that orbit drag can s ign i f icant 1 y affect the solar synchronization 
and eclipsing history of the satellite's orbit. Another consideration which 
will be discussed in the next section is the effect that orbit decay can 
have on the earth coverage pattern , i.e., the way in which the orbit tracks 
are phased so as to cover a desired mapping area on the earthd surface. 

As will be seen, this is probably the most critical aspect when considering 
how much orbit decay, if any, is acceptable. 

Orbit drag can be virtually eliminated as a consideration over a typical 
satellite's lifetime by choosing an altitude of perhaps 800 km or more. This 
may be preferred to avoid the cost, complexity and weight associated with 
having to provide a velocity compensation system for drag makeup during the 
lifetime of the satellite. On the other hand, if satellite RF power is a 
critical factor it may be desirable to go to a lower orbit to conserve power 
in which case the satellite could either be allowed to decay with no drag 
compensation or a velocity correction system could be provided. Allowing 
the satell i te to decay with no drag compensation would be acceptable only 
if the amount of decay over the satellited lifetime was not in any way 
detrimental to its operation and mission. As will be seen in the next 
section, the requirement to maintain a pre-determi ned earth coverage pattern 



21 


would allow virtually no orbit decay over the sate 11 i te*s lifetime. The 
question then becomes one of determining how low an orbit might be 
feasible assuming a drag compensation system is included on the satellite. 

Orbit sustaining could be accomplished in two ways. The first method 
is simply a continuous application of thrust equal to drag. The second 
method consists of altitude gaining Hohmann transfers in between periods of 
altitude decay due to drag. Thrust is applied for only short periods of 
time and impulsive applications are assumed. In this method the altitude 
is allowed to vary an amount Ah whose upper limit is the desired altitude 
The continuous thrust case would require an electric propulsion 
system such as an ion thruster or resistance jet. It can be shown 
that the velocity requirements for continuous thrust orbit 
sustaining are given by 



where 


I = the specific impulse (sec) of the rocket 

SP 2 
g = the acceleration of gravity, 9.80^65 m/sec 

= time (sec) divided by the ballistic coefficient (defined previously) 


W/C q A 


2 

q = the dynamic pressure, kgf/m 

A plot of equation (9) is shown in figure 10 for various initial 
altitudes and three representative specific impulses. Note that for altitudes 


of 500 km or higher the curves for the l gp range of 1000-3000 sec are very 
close so that only the 1000 sec curve is plotted for clarity. 



22 


The thrust required to cancel drag as a function of altitude is shown 
in figure 11 along with the propellant consumed over a 5 -year period for 
l s p>*s '000, 2000 and 3000 sec. Note that the required thrust ranges 
from about the size tested on SERT II down to the microthruster range. 

Two cesium contact ion microthrusters with several discrete thrust 
levels ranging from 4-20 A' b were tested on ATS IV in 1968 (ref. 7 ) and 
generally performed as expected. 

It can be seen from figure 11 that the amount of propellant consumed 
over a 5 -year lifetime would be reasonably small for altitudes of about 
400 km or more. At lower altitudes the propellant could represent a sig- 
nificant fraction of the satellite’s weight which Is undesirable since It 
subtracts from useful payload. The feasibility of using an electric pro- 
pulsion system for drag compensation and its impact on the satellite's design 
cost, complexity, reliability, etc. requires much more detailed study than 
is possible in this report. The major problem would seem to be one of 
developing a sufficiently high level of confidence in the ability of the 
system to operate effectively over the satellite's lifetime to warrant its 
use . 

The impulsive thrust case would require a chemical propulsion system 
such as a monopropellant hydrazine system. The velocity requirements in this 
case are determined by the Hohmann transfer maneuvers from the low altitude 
to which the satellite has been allowed to decay to the original starting alt 
tude. Two impulsive velocity corrections are required for each maneuver : 
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(1 ) A V j at the low alti tude to place the satellite on a transfer ell i pse with 
an apogee at the original desired altitude and (2) AV^ at the apogee of the 
transfer ellipse to recircularize the orbit at the desired altitude. The 
required velocity increments are given by 



where 


(10a) 


(10b) 


r. = radius of initial orbit at beginning of Hohmann transfer maneuver 

r^. = radius of final desired orbit 

The Hohmann transfer maneuvers are separated by periods during which 
the orbit is allowed to decay to some lower altitude limit so that the 

thrusters operate on a cyclic basis. The duty cycle depends on the starting 

altitude (drag encountered) and on the amount of decay allowed. 

A plot of the total AV (AV. + AV ) versus rj/T — is 

1 Z W/ uA 


presented in figure 12(a) for various initial altitudes and for A h/h Q = 

.10, .05 and .01. and Av^ are near ly equal due to the proximity 

of the upper altitude limit (desired orbit) and the lower limit (decayed 

orbit). The duty cycle structure of these curves was omitted for clarity 

since the number of cycles gets fairly large for low altitudes. A plot of 

the duty cycle or time between Hohmann transfer maneuvers is shown in figure 

12(b). As an example, for ho = 400 km and A h/h Q - .05 the duty cycle would 

2 

be t = 0.53 days so that for a ballistic coefficient of 50 kgf/m 
W/C^A kgf/m 2 
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the Hohmann transfer maneuver would be made every 26.5 days. Each maneuver 

involves two thruster firings so that for a 5-year mission the thrusters 
would fire about 140 times. For higher altitudes the duty cycle is increased 

and the number of thruster firings would be reduced as indicated in figure 
12(b). 

In comparing figure 12(a) with figure 10 it may be seen that the 
AV's for either system compare reasonably well. However, since the 
chemical system has a much lower specific impulse, the propellant require- 
ment will be greater (compare figure 13 with figure 11). 

If we were to limit the propellant consumed for drag compensation to 
say 5 percent of the initial satellite weight for a 5-year mission the minimum 
initial altitude allowable would be about 525 km for a 300 sec thruster and 

Ah/h = .05 (see figure 13)* 
o 

3. Earth Coverage Patterns 

The considerations discussed thus far for sun-synchronous orbits have 
typically been in the nature of constraints bearing upon the choice of an 
orbit. These constraints help to set a lower limit on the choice of 
altitude but are not really definitive. Theoretically, an infinite number of 
orbit altitudes exist from which to choose so that it is desirable to have 
a more definite criteria upon which a reasonable choice can be made. For 
a surveillance and mapping satellite one of the important criteria involves 
the way in which the ground areas of interest are covered or scanned. This 
involves consideration of the distance between successive or neighboring 
•^not necessarily successive) orbit tracks across the surface of the earth, 
and the time required to completely cover an area. 
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To focus this discussion, the application of surveillance and mapping 
of ice on the Great Lakes will be considered„ The Great Lakes are bounded 
roughly by latitudes of 41 deg north and 49 deg north and longitudes of 75 
deg west and 92 deg west. This area is shown in figure 14 and would measure 
about 1400 km (southern edge) by 1200 km (northern edge) by 900 km (meridian). 

The orbit will pass over this area twice a day, once on the up leg 
(south to north on a north-northwesterly heading) and once on the down leg 
(north to south on a south-southwesterly heading). However, only one of 
these passes (either one) can be considered for an orderly mapping since the 
two passes will cross the area out of sequence and at different angles. 

Only the up leg (south to north) will be considered here. 

Once the orbit is established, the orbit track on the earth‘s surface 
will regress or move westward as the earth rotates. For the limited range 
of orbit altitudes being considered, the equator crossing will regress about 
.38-. 45 rad/orbit (22-26 deg/orbit). For example, if we had an orbit with 
a period of 1.5 hours (altitude = 274,400 km) the equator crossing would 
regress 0.3927 rad/orbit (22.5 deg/orbit) since the earth rotates 0.2618 
rad/ hour (15 deg/hour). In this case, exactly 16 orbits would be completed 
in one day with the 1 7 th orbit repeating the 1st. 

Obviously, if we reduce the number of exact orbits completed in one day 
to 15* 14, etc, the orbit period will become larger and the amount each 
successive track moves westward will increase. In each case the required 
orbit altitude is unique. For example, 14 orbits/day requiresan orbit period 
of 1.7143 hours and an altitude of 893*777 km. The orbit would regress 
0.4488 rad/orbit (25.714 deg/orbit). 
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It is apparent that the separation between successive orbits (orbits 
1 and 2, 2 and 3» etc.) would be at least 0.3927 rad (22.5 deg) In longitude. 
Thus as successive orbit tracks pass over the Great Lakes, for example, they 
will be separated by at least 1650 km minimum which is wider than the entire 
area. The radar would thus be required to cover a swath of 1650 km minimum 
to obtain an image once a day. Initial study indicated this was not feasible 
for a reasonably sized multi-beam antenna. This is discussed in more detail 
in the radar section to follow. 

Smaller orbit track separations can be obtained by considering patterns 
that repeat every 2 days, every 3 days, etc. The longer the cycle, the 
smaller the track separations will be. The pattern need not repeat in an 
integral number of days so that there are virtually an infinite number of 
repetition patterns with their corresponding altitudes to choose from. 

However, it Is not necessary to consider all possibilities since consideration 
of patterns that repeat in an integral number of days yields a reasonably 
large number of possible orbits from which to choose and simplifies the analys 
cons i derably. 

Since the earth rotates 2 A 1 radians in one day the nodal rate of change 

due to the earth's rotation can be found from 
¥ ^ 

fl = 2 'ft ± f b . 

J * a a rad/orbit 

D 

N 

or A = 2 7f D rad/orbit (11) 

(DN ± 1 
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where D =' number of days between pattern repetitions 

N = approximate number of integer orbits/day (rounded to next 
highest integer orbits/day) 

Note that although an integer is used for N to generate a set of 

A _ the number of orbits actually completed in one day will not 
a 

necessarily be an integer . Also, the plus sign in the formula gives a 
nodal precession (eastward movement) while the minus sign gives a nodal 
regression (westward movement). 

As an example, suppose we want the orbit track pattern to repeat in 
a westerly direction every two days (D — 2) with N = 15 orbits/day. Then, 


A = 47f = 0.^3332 rad/orbit 

A a 29 

Since the earth rotates 2 IT rad in 24 hours the required orbit period 
would be 0.43332 X 24 = 1.6552 hours. The required altitude would be 

W 


725.630 km. For these conditions 14.5 orbits would actually be completed 
in 24 hours with 29 orbits being completed in 48 hours. The 30th orbit 
would exactly repeat the track of the 1st orbit and the cycle would be 
repeated. The separation distance between the 1st orbit and the 15th 
or between the 2nd and 16th and so on would be 0.21666 rad in longitude 
or about half of what it would be for 15 orbits/day repeated daily. 

It should be noted here that since the orbits under cons iderat ion 
are required to be sun-synchronous they would have a small precession rate 
about the earth due to earth's oblateness. However, this precession rate 
does not affect the preceding calculations since it is accounted for 
by the motion of the earth about the sun. The length of the day is 
accordingly reckoned at 24 hours to account for this motion. 
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A partial map of sun-synchronous orbits is shown in figure 15 where 
orbit track separation in longitude is plotted versus altitude. Also 
indicated on the figure for some discrete points are the exact number of 
days before the orbit track pattern repeats itself and the number of the 

orbits completed in that time period. 

In the case of surveillance and mapping of ice on the Great Lakes 
where the conditions being surveyed can be quite changeable due to current, 
winds, temperature, etc. it would be necessary to obtain a complete image 
of the region in as short a time as possible to effectively aid navigation 
and ice control operations. In fact short response time may be necessary 
for many applications (ocean state surveillance, pollution monitoring, etc.). 
Short response time infers wide swath width coverage. As mentioned earlier, 
however, the radar system would be limited in the maximum swath width 
achievable which at least rules out obtaining a complete image daily of an 
area the size of the Great Lakes from a sun-synchronous orbit. For example, 
if we consider an orbit track separation of .1 rad in longitude, the orbit 
track separation or swath width at the latitude of the Great Lakes would be 
about 450 km. If this were the maximum achievable swath width, from figure 
15 this would limit the choice of orbits to those having a track repetition 
cycle of about 4 days or greater. However, the track separation distance 
at the equator will be about 45 percent greater than at a latitude of 45 
deg. Thus, if the satellite were required to map any part of the globe, 
the equatorial separation distance would be the design point. If the 
maximum separation distance were limited to 450 km at the equator it would 
take about 6 days or more to obtain a complete image of the Great Lakes. 
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The separation between orbit tracks is a function of the orbit 
inclination, the longitudinal separation and the latitude being con- 
sidered. The separation distance will be maximum near the equator and 
zero at the point at which the two orbit paths cross which will tend to 
be near the maximum latitude point of the orbits □ In order to compute 
the separation distance as a function of latitude it is necessary to 
choose a specific sun-synchronous orbit. A method for calculating 
the separation distance between any two orbit tracks on the earth's 
surface displaced in longitude by an amount A is presented in 
Appendix B. 

C. Equatorial Orbits 

This class of orbits is of interest since it allows the earth scan 
velocity to be small by choosing an altitude somewhat below geosynchronous 
so that the satellite has some desired easterly velocity relative to the 
earth's surface. As will be discussed in the radar section this 
would enable the radar system to obtain high image 
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resolution while covering a large swath width with only a few or even a 
single beam. This type of orbit would experience no atmospheric drag which 
is a simplifying feature. Another advantage than an equatorial orbit would 
have over a sun-synchronous orbit is that areas can be scanned more than 
once per day. On the other hand, orbit eclipsing will be more prevalent 
than with the sun-synchronous orbits. Also, limitation on the radar beam 
incidence angle will preclude mapping the near equatorial region (angle 
too vertical) and the polar regions (angle too shallow) which may be a draw- 
back for some applications. It may be possible to eliminate this problem 
by giving the orbit a slight inclination to allow the side-looking radar 
to periodically view the equatorial region while not significantly increasing 
the satellite's velocity with respect to the earth's surface. This was not 
investigated in this study but should be considered in any future, more 
detailed studies of this concept. 

1 . Orbi t Eel i pse 

The geometry of a circular equatorial orbit with respect to illumination 
from the sun can be considered in basically the same way as for the sun- 
synchronous orbits discussed earlier. In this case, however, the orbit 
inclination is zero and there is no sun-synchronization. Maximum daily 
eclipsing will occur when the sun is crossing the equator on about March 21 
and again six months later. Shown in figure 16 are plots of how the eclipse 
time varies throughout the year for a range of equatorial orbit altitudes 
from 10,000 km to 40,000 km. Geosynchronous altitude is 35*786.04 km. 

2„ Earth Coverage Patterns 

Earth coverage from an equatorial orbit depends primarily on the 
satellite's velocity with respect to the earth's surface and the latitude 
being scanned. The relative angular velocity, a? r > between the satellite 
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and the earth is 

> 

Ui = U> ± Lu> 

r s e 

where is the satellite's angular velocity 

is the earth's angular velocity 

The relative velocity of an equatorial satellite with respect to a point 
on the earth's surface is 

V = U X R = u; R cos V (12) 

rre re 

and is plotted versus altitude in figure 17 for various latitudes. For 
altitudes below geosynchronous, the difference between the satellite's and 
the earth's rotation will cause the satellite to "lap" the earth so that 
multiple daily scans may be achieved. Shown in figure 18 is the number of 
daily passes over an area on the earth's surface versus satellite altitude. 
It can be seen that areas to be scanned may be covered at least daily for 
altitudes below about 20,000 km (1/day occurs at 20 , 207.8 km). 

The number of satellite eclipses encountered each day, N , will occur 
once more than the number of target passes per day. The longitudinal 
separation of these eclipses (center to center) would simply be 360/N^ 
deg. The longitudinal slice "blacked out" on each eclipse measured as a 
segment of earth's longitude would be 
360 f (24- P) deg 

where P is the satellite's period (hours) 

f is the fraction of the orbit in shadow (figure 16) 
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Tor example, for a 4/day coverage where h = 8062.8 km, P - 4.797** 
hours and f = .1465 (max.) each eclipse would cover 42.2 deg of earth's 
longitude with each eclipse being separated (center to center) by 360/5 = 

72 deg. The longitudinal location of the eclipses would depend on the 
initial relationship between the time of day and the longitudinal position 
of the satellite with respect to the earth. For an integral number of 
passes per day the eclipses would always occur at the same longitudinal 
positions. [n the above example all longitudinal positions would see full 
sunlight at least three times per day. 

For a satellite in equatorial orbit, the radar beam incident angle 
will depend on the orbit altitude and the latitude being mapped. This is 
illustrated in figure 19. If radar system performance were limited to 
incident angles of about .35 rad (20 deg) for instance, latitudes above 
60 deg could not be mapped. However, even this amount of coverage would 
cover most of the world's shipping lanes, for example, suggesting an 
application for ocean state surveillance and ice imaging. 

Of importance to the radar system is the rate at which the distance 
to a target varies with time in order to obtain a Doppler effect. This 
distance r ( (J) ) is shown in figure 20 and r ((J) ) may be found as follows; 

From figure 20, beam position 1 lies in the plane of the subsatellite 
meridian at some latitude and represents the point at which r ( $ ) is 
minimum. Point 2 represents a beam position ahead of the satellite. As the 
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satellite moves in an easterly direction the radius from the satellite 
to position 2 will become shorter until it reaches the value represented 
by position 1 and then ft would lengthen again as the satellite moves 
further eastward. 

From the law of cosines 
2 2 2 

r = R + R - 2 R R cos £ 
s e s e 

where cos £ = cos V cos 

Differentiating with respect to time we get 
2 rr = 2 R s cos X sin 

but (p = U> 

Then r = R g R^ to cos X sin ^ (13) 

r 

22 , ] /2 
where r=(R +R -2R R cos * cos $ ) 
x s e s e ~ 

3. Radiation Damage to Solar Cells 

The existence of trapped, high energy electrons and protons and solar 

flare protons in the equatorial region about the earth can cause significant 

damage to unprotected solar cells and other electronic components. Electronic 

components can generally be shielded to protect them from radiation damage 

with some weight penalty. However, solar cells present a more difficult 

problem. 

Calculations were made to determine the altitude and cover 
glass thickness necessary to achieve a five year lifetime with no more than 
a 25 percent degradation in solar array power. A nominal 10 Ct -cm n-on-p 
cell with a thickness of 12 mils and infinite back shielding was assumed. 

It was determined that an equivalent 1 MeV electron fluence of approximately 
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14 2 

7 X 10 e /cm would reduce the cell power by 25 percent. Accordingly, 

for a five year mission the maximum annual fluence impinging on the cell 

14 2 

that could be tolerated would be 1 .4 X 10 C /cm . 

Equatorial orbit altitudes ranging from 1000 to 20,000 km were studied. 
For orbits below 11,000 km the trapped protons were the dominant damaging 
particles, while trapped electrons were predominant above 11,000 km. 

Solar flare protons made a significant contribution to the solar cell 
damage only at altitudes above 14,000 km. The solar flare model with the 
greatest proton fluence was chosen. It was concluded that a five year life- 
time can be obtained only for orbital altitudes below 1200 km or above 
11,000 km. The heart of the trapped proton belt lies between these altitudes 
and lifetimes less than a year would be expected. Outside these extremes, 
the effect of cover glass thickness was considered and it was concluded that 
a five year lifetime with no more than a 25 percent power degradation could 
be achieved with cover glass thicknesses less than 30 mils in all cases. 

However, equatorial orbits below 1200 km are not practical since earth 
coverage would be severely limited by radar beam incident angle constraints 
at low altitudes. For example, at 1200 km, earth coverage would be limited 
to a narrow swath between 10-20 deg latitude. The region of interest for 
a solar array system would thus be confined to altitudes above 11,000 km. 

The altitude should be chosen no higher than necessary since the required 
RF power increases fairly rapidly with altitude. This will be discussed in 
more detail in the following section on Radar Considerations. 
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Within the 1200 to 11,000 km region another type of power system 
could conceivably be used such as an RTG (Rad io- I sotope Generator). Other 
electronic components would need special radiation protection thus incurring 
some weight penalty. Such a system is judged to be more complex and 
costly than a solar array system operating outside this region. However, 
this approach would need more detailed study than can be given here to 
determine its feasibility. 
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RADAR CONSIDERATIONS 

Ag General 

Preliminary results of the Lewis airborne SLR ice mapping program 
indicate that rapid coverage of the Great Lakes is necessary (ref. 2), 

Figure 21 shows typical imagery of Lake Erie obtained with that system. 

The figure shows a complete pass on February 22, 1973. The dark line 
through the center of the image is the blind ground track of the aircraft. 

A portion of the February 23, 1973 pass is also shown with a sketch of 
Lake Erie. Observe the February 23 pass and note the ice pack against the 
western edge of Pelee Point. On the 22nd this was relatively open water. 
Close scrutiny of these photographs will reveal many other striking changes. 

This same phenomenon has been observed with ERTS-1 imagery an example 
of which is shown in figure 22. Note the change along the southern shore 
that occurred in one day. 

Because the Great Lakes ice formation is a dynamic phenomenon, rapid 
coverage of the lakes is desired. Daily coverage would be desirable and 
more than 3**4 days separation would render the imagery useless if accurate 
navigational maps are required. However, these lakes span approximately 
8° in latitude by 17° in longitude. This is an area of about 900 km by 
1300 km. Rapid coverage of such a vast area will require very large swath 
widths. As a result of antenna dimensional limitations and signature 
constraints, spaceborne synthetic aperture radar (SAR) proposals to date 
limit expected swath widths to about 60-100 km. Rapid coverage of these 
lakes will require swath widths 4-10 times greater. Among the obvious ways 
of obtaining larger swath widths are: adjustment of the SAR look angle 

(angle from vertical as seen by the spacecraft) and/or using larger antennas. 
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Since target signatures are sensitive to look angle, a look angle 
favorable to large swath widths may not be favorable to target signature. 
Therefore, it is desirable to have flexibility in choosing the look angle. 

Since the realizable swath width is proportional to the antenna 
azimuth dimension, using larger antennas is a conceptually simple alternative. 
Deployable antennas of 10 meters in dimension exist (ATS-F satellite) and 
larger antennas may be possible. A 10 meter antenna would allow a swath 
width of about 60-100 km at low altitudes ( 1000 km) and for small 
incidence angles 10 o -40°). A 4-10 fold increase would imply a deployable 
antenna of 40-100 meters. Such a development for X-band (most favorable 
band for ice imaging) would be a formidable task. Hence it would be 
desirable to find other schemes. 

A possible solution has been proposed (ref. 8) which will yield large 
swath widths with reasonable antenna size and no restrictions on look angle. 
Essentially the method employs independent SAR's in parallel to increase 
the realizable swath width. This section will present a simplified analysis 
of this scheme and use it to estimate spacecraft requirements for wide swath 

i 

coverage . 

B. Synthetic Aperture Radar 
lo Single Beam Radar 

The typical spaceborne SAR geometry appears as shown in figure 23. As 
the spacecraft moves along its trajectory, RF bursts are periodically 
directed at the ground. The period of the bursts and the spacecraft velocity 
relative to the ground target determine the separation, <!iX, between 
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transmit/receive locations as seen by the ground target. The beam 
strikes the ground to one side of the satellite radar and a portion of 
the reflection returns in the direction of the spacecraft. This side 
looking geometry is necessary for obtaining good range resolution since 
the targets on the ground can be discriminated according to time delay. 

For example, the side view of figure 24 shows that the reflection 
along lags the reflection along Rj . (f the transmitted pulse width is 
T then the minimum resolvable target separation, A r, 

A'f > c r 

Z ( - -©■ 

where c is the velocity of propagation and 0 is the local incidence angle. 

Therefore, resolution in range is a matter of choosing an appropriate pulse 

width and adjusting the transmitter and receiver bandwidths accordingly. 

The transmit/receive locations along the trajectory offer different 

radar "views" of each target. If the returns at each point are recorded 

and later added together vectorial ly, the result is equivalent to the output 

of an array of antennas with the same spacing and total length. The apparent 

or synthetic length parallel to the flight path, d , and the dimension 

as 

perpendicular to the flight path,. D , form the synthetic aperture. The 
null to null azimuthal resolution of this SAB, Aa, is 

cc > A 

where \ is the wavelength of the incident radiation. The equality is 
obtained for uniform weighting of the radar pulses and perfect coherency 
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between successive pulses. Immediately, one can see a factor of two 
improvement over conventional arrays. This is an advantage obtained 
through coherent detection. However, other 1 imi tat ions ar i se which make 
the equality unattainable in certain applications. Cutrona (ref. 9) has 
shown that the curvature of the reflected waves produce phase errors which, 
when unaccounted for, limit the resolution. If unaccounted for, these 
phase errors limit the resolution to 



2 

where R is the slant range to the target. However, if compensation is 
made for these phase errors the resolution is theoretically (ref. 9) 

A®- - 

2 

where D a is the real antenna azimuth dimension. However, perfect coherency 
between pulses is not always possible and atmospheric variations become 
significant when d is long (ref. 10). Nevertheless, resolutions of 

3 S 

10's of meters are obtainable from altitudes of 1000's of km. 

The performance of such antennas has been analyzed thoroughly and 
results can be found in the literature (refs. 11, 12, and 13). The gain 
of such an array is given by 

G = N 2 G 

o 

where N is the number of RF bursts coherently summed and G q is the real 
antenna gain. For non-uniform weighting some loss in gain can be expected 
but sidelobe suppression is achieved (ref. 12). 
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Conventional linear arrays have only a slight effect on the effective 

receiver noise temperature. SAR's do not have this advantage. The gain 

of the SAR is the result of adding many radar returns. Unfortunately, the 

noise added to each return by the receiver front-end will contribute to 

each summation. Consequently, as successive range returns are coherently 

added the signal power goes up as the square of the number of pulses and 

the noise power goes up proportional to the number of pulses. Therefore, 

the net improvement in signal to noise ratio only goes up proportional to 

the number of pulses. This effect will be accounted for when estimating 

the average SAR power requirements. 

Since the generation of the synthetic aperture is periodic, major 

sidelobes occur at intervals of . radians, where 

maj 

Sl . = A- 

maj Z AX 

These major sidelobes must be suppressed by the real antenna since the 
synthetic weighting, N , is the same at all major sidelobes. This places 
a restriction on the minimum size of the real antenna or maximum element 
spacing through the inequality, 

A. > x. 

ZAX " Dc. 

where 0 is the real antenna azimuth dimension. Clearly, the constraint is 
a 

independent of wavelength and can be interpreted as 
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A plot of a typical receiver pattern is shown in figure 25. The dotted 
line indicates the response of the real antenna and the solid line is the 
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overall response. The location of the real antenna nulls and spurious 
lobes are indicated. Clearly, by adjusting the real antenna azimuth 
dimension or the synthetic array spacing, the major sidelobe can be moved 
to a minor response of the real antenna. 

Therefore, by choosing an appropriate pulse width, *£ , and synthetic 
aperture, d , a resolution cell 

3 5 

AS 3 Ah* a. a- = JBA £* 

A das CiS -Gi- 
ls achieved (provided > Ql* for a focused SAR and x \Ja£ 

'2. TTs ~T 

for an unfocused SAR). 

Processing techniques exist for generating the vector sum of N returns 
at all range intervals simultaneously. The favored technique presently is 
an optical process (refs. 14, 15, and 16) which reduces each of N range 
scans to a single scan of high resolution. The high resolution scans are 
then recorded side by side on photographic film to produce a high resolution 
radar image of the ground. 

To avoid range and azimuth amb i gu i t ieSj constra i nts must be placed on 
the radar system parameters. The azimuth ambiguity problem has already been 
mentioned. The requirement was finally interpreted as a restriction on 
element spacing 

/A* Q&. 

2 

Since the element spacing is fixed by the radar pulse repetition frequency, 
prf, and spacecraft relative velocity, V , the above inequality can also be 

P rf 2 2 v 

r 
D 

a 


written 
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Therefore the minimum prf is established to avoid azimuth ambiguities. 

2. Single Beam/Multiple Frequency SAR 

As was mentioned before it is desirable to have as wide a swath width 
as possible. The swath width is fixed by the prf. This can be seen by 
examining figure 26. The reflection along R 2 lags the reflection along 1^ . 
If the prf is high enough the reflection along R 2 for pulse p will interfere 
with the reflection of pulse p + 1 along Rj . This causes a range ambiguity 
which limits the prf to 

prf ^ 

Z (R z " Ri) 

However, consider the case when successive pulses are transmitted at 
different frequencies. If these frequencies are spaced so that no inter- 
ference occurs between carriers, then the above restriction only applies to 
eadjcarr ier individually. But if the overall pulse rate is prf then the 
pulse rate for each carrier is prf/m where m is the number of carriers. 
Therefore, the overall prf limitation for mul t i frequency operation becomes 


prf _ 

The prf is then constrained by range and azimuth ambiguities to lie between 



0 ^ 4 " 

;4CKx-Ki) 


a 

3. Mul t i frequency/Mu 1 t i beam SAR 

Range discrimination can also be obtained by using several focused beams. 
The above restriction on prf would then apply to each beam. Sidelobe inter- 
ference could be suppressed by using different bands of frequencies for 


each beam. 
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It is clear from the above equation that the widest possible swath 


occurs for 


prf = 2 V r 

D~" 


a 

The largest possible unambiguous range difference is then 

R 0 - R, = m C 
Z 1 prf 

Referring back to figure 26 one can see the geometrical relationship 
between beamwidth, /3 , swath width, SW, and range difference, 4 R. For 
a given angle <fi between nadir and the target the slant range is given by 



where R.tradius of the earth 
e 

R^ spacecraft orbital radius 

4 > is chosen to produce a favorable incidence angle, 8j . 
particular incidence angle, 0^, ^ ^ can be computed from 
P 1 - cos" 1 / R e cos 0 ] 

V s 

Letting ^ be the smallest angle from nadir, then 



Given a 


R 


1 
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The range difference is obtained from the equation for prf 




and therefore, ft, - 4 R 




The most extreme angle, from nadir is then computed from 

(t = r •-'$_[_/ .4- Ke _ r-J- \? 

*2 Co.- (a l -Ci?e)3 

The resultant swath width is then 

S' ^ ~ K e C ^ - (4, ) 


and the required antenna range beamwidth is 



Therefore, the required antenna range dimension is 

n _ 2 a 

' r j3 


For multiple beam operation the angular width of each beam is fixed 
by the antenna range dimension. The side looking geometry causes beams 
farther from nadir to cut wider swaths than those nearer nadir. Therefore, 
the range ambiguity problem is most severe at extreme ranges from nadir. 
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If D ^ is chosen to suppress the range ambiguity at ranges distant 
from nadir, then the beam close in to nadir will be too narrow to realize 
the full unambiguous range interval. This effect can be compensated for 
by increasing the number of independent carriers in those beams extreme 
from nadir. 

Letting SWj be the swath width of the j th beam the total swath width 


is then. 


s-w- Z S ICj 
j 


where b is the number of beams. 

Assuming processing in all 

length, d , is the same for all 
a s 

beam becomes 

A o-j - 


beams i synchronized so that the synthetic 
beams, the azimuth resolution of the j th 



Choosing d as to give a nominal resolution in some beam between j = l and 
j * b, say j - n, then the azimuth resolution in the j th beam can be 


wr i tten 




A c\ 


where A a is the nominal azimuth resolution and is the slant range of 
the n th beam. 

Assuming the radar pulse width is the same for all beams, the range 
resolution in the j th beam is 

A V‘ - — A-,,.. 1 : 
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where 6 ^ is the local angle of incidence of the jth beam (approximately 
the same across the beam for swath widths encountered in this study). 
Again, choosing / t to give a nominal range resolution, Af , in the nth 
beam, the range resolution in the jth beam is then 

A - A V Cor.. Gn_ 

6^ 


In the jth beam the received signal power can be written. 

Pc 


where 
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is the received signal power in the jth beam 
is the RF losses in the jth beam 

is the transmit antenna gain in the jth beam direction 
is the space loss in the jth beam 

is the scattering coefficient in the jth beam direction 
is the effective gain in the receive mode 
This power must exceed the receiver noise power by the desired signal to 
noise ratio. Therefore, 
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where KTB is the effective receiver noise power 

N is equal to the number of RF pulses used to generate the 
synthetic aperture 

(S/N) is the desired signal to noise ratio 
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where P tJ c average transmitter power in jth beam 

t - efficiency of real aperture 
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the total required average power can be written 

p t =. )' 4ii Ky ft," W- i i C K r c C± 

D«f D-,’ V/ Ay A <r c os &n 




where t> , the number of beams satisfies. 
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The RF bandwidth of the radar may vary from beatn to beam since some 
beams may have multiple carriers. The video bandwidth is assumed to be the 
saime for each beam and this was used in estimating the average power require- 
ments. However, the radar front end and the data link must pass all beam 
returns simultaneously. Therefore the RF bandwidth of each radar beam is 

- Z 

'T 

where m. is the number of carriers in the beam and X is a factor 
J 

accounting for any guard bands used. The data link requirements are not 
included in this study. 

C. Application of Mul t i f requency/Mul t i beam SAR to Equatorial and Polar 
Orbits 

The relatively small swath widths of conventional orbital SAR's is a 
result of the high velocity required to orbit. When a target comes into 
the field of view of the real beam pattern, it produces a doppler shifted 
return depending on the length of the synthetic aperture and the relative 
approach velocity toward the target. The faster the spacecraft the higher 
the frequency. The sampling theorem dictates that sampling of the doppler 
frequency must exceed twice the doppler frequency. Thus the prf increases 
as the spacecraft velocity increases. A high prf reduces the obtainable 
swath width through the range ambiguity. Therefore, as the velocity 
increases the obtainable swath width decreases. This property is displayed 
in figure 27. The incidence angle limit contour represents the distance 
from nadir to a 10° incidence angle. The swath width for lower velocities 
encounter this limit. However, for orbital velocity the obtainable swath 
width is far below this limit. 


49 


1 . Equator i al Orbi ts 

Orbital velocity decreases as the spacecrat altitude increases. In 
equatorial orbits the angular velocity of the earth subtracts from the 
angular velocity of the satellite creating an apparent low orbital velocity. 
Hence, more improvement in swath width should occur for equatorial orbits 
than for polar orbits. 

Unfortunately, striking improvements are not obtained until very high 
altitudes are used as shown in figure 28. These contours were generated 
for a target at 48° latitude (northern extreme of the Great Lakes) and an 
X-band SAR. Two azimuth dimensions were used. For a 10 meter azimuth 
dimension the solid line results. At each altitude the antenna range 
dimension is set equal to that required to eliminate the range ambiguity. 

For altitudes below approximately 11,500 km, a 10 meter range dimension 
is insufficient to eliminate the range ambiguity. Two frequency operation 
is required to enlarge the range interval to & value controllable with a 
10 meter range dimension. We are interested in the maximum obtainable swath 
width so the range dimension is reduced to the minimum necessary to eliminate 
the range ambiguity. Above 11,500 km, the 10 meter range dimension limita- 
tion is sufficient to eliminate the range ambiguity and mul t i frequency 
operation is not required. For a 15 meter azimuth dimension, single 
frequency operation is realized over the full altitude range. The limiting 
contour on the left side of the figure is the 10° incidence angle limit. 

The limiting contour at the top of the figure is an arbitrary 900 km limit 
built into the design procedure. 900 km is sufficient to span the five 
Great Lakes in a north-south direction and a larger swath is not necessary. 
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At these high altitudes severe requirements are placed on the space- 
craft power supply. The required average power is shown in figure 29. 

For a 10 meter azimuth dimension, the required range dimension for multi- 
frequency operation below 11,500 km is small and causes a large increase 
In required power. Above 11,500 km, there is a drop in the unambiguous 
range interval and a corresponding increase in the antenna range dimension. 
Therefore, the required power drops. With a 15 meter azimuth dimension 
the required power is monotonic to about 24,000 km. Between 11,500 km and 
approximately 24,000 km there is very little difference in the required 
power for the two antennas. However, as was shown in figure 28, the swath 
width is 50% larger for the 15 meter antenna. 

At altitudes near synchronous, a M hole M in the power requirement 
occurs. In generating the data of figures 28 and 29 an arbitrary limit 
of 900 km was placed on the swath width. At altitudes above 28,000 km, 
the single frequency range ambiguity exceeds this limit for either antenna. 
Therefore, the prf was increased to fix the range ambiguity at the 900 km 
limit. Normally the prf would be chosen as small as possible to maximize 
swath width. In this case the swath width is limited by another constraint. 
The extra returns can then be used to increase the SAR gain and less power 
is required. This reduction is most noticeable near synchronous where the 
relative velocity of the spacecraft approaches very low velocities. The 
required synthetic aperture is fixed by the resolution requirement and remains 
approximately constant. The prf remains fixed to specify the 900 km range 
ambiguity. Because of the decreasing velocity more and more excess returns 
are available. Consequently, the resultant increase in SAR gain causes 
the power requirement to experience a rapid drop near synchronous altitude. 
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As synchronous altitude is approached, the orbiting time of the 

spacecraft becomes very long and the rate of coverage so low as to make 

such a scheme impractical for surveillance of the Great Lakes. Another 

serious weakness of these orbits is the required stability of the SAR 

local oscillator. Figure 30 depicts the rapid rise in the required long 

term stability as the spacecraft approaches synchronous orbit. Stability 

1 2 

in excess of 1 part in 10 places a very severe restriction on the local 
oscillator. 

An examination of figures 28 and 29 reveal that a single pass coverage 
of the Great Lakes (900 km swath) from equatorial orbit has a very demanding 
power requirement. However, it is interesting to note that for equatorial 
orbits every pass is a potential imaging pass. Thus, by constructing a 
mosaic of .five consecutive passes, for example, a daily map of the Great 
Lakes ice cover is obtainable with an X-band SAR at an altitude of approx- 
imately 6400 km. If a 15 meter antenna is used, single frequency operation 
is sufficient with an expenditure of about 370 watts per pass. For a 10 
meter antenna, two frequency operation is necessary with an expenditure of 
1.2 kw per pass. Clearly, development of a 15 meter X-band antenna would 
be a tremendous asset to such a system. In either case, the local oscillator 
stability requirements are modest and approximately 2 parts per 10^. It 
should be noted that at 6400 km altitude, latitudes to 50 ° are accessible 
witha 10° incidence angle (figure 19 ). 

The question arises as to the applicability of such a system at lower 
latitudes. At lower latitudes the incidence angle is greater and for 100 
meter resolution, the required pulse length must be smaller. This requires 
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a corresponding increase in receiver bandwidth and hence, transmitted 
power. Simultaneously, the achievable swath width increases requiring a 
corresponding decrease in the antenna range dimension. This reduces the 
gain of the antenna and the required power is further increased. This 
trend is depicted in figure 31 . 

If the achievable swath width at the lower latitudes is greater than 
required, a higher prf and larger antenna range dimension can be used to 
reduce the required power. For example, the achievable swath width at 10° 
latitude for an X-band SAR at 6400 km altitude (figure 31) is 409 km. The 
required power is 7 kw. Quadrupling the prf reduces the swath width to 
100 km. The antenna range dimension goes from 1,06 meters to 4.24 meters. 
The required power is reduced by a factor of 16 to 426 watts. 

Both the Great Lakes application and low power low latitude applica- 
tions could conceivably be satisfied by using an antenna with multiple 
beams. That is, a large number of beams could be used at the lower 
latitudes with a smaller number (down to a single beam) used at the higher 
latitudes. This would be equivalent to having an antenna with adjustable 
range coverage. 

The multi beam antenna approach would allow tailoring of the SAR 
characteristics to optimize power and performance at each latitude. Hence, 
the equatorial SAR would be applicable to the lower latitudes as well as 
the latitude range of the Great Lakes. 

While an altitude of 6400 km for an equatorial orbit appears to be a 
reasonable compromise between the power required, achievable swath width, 
antenna size and frequency of coverage, the satellite would be in the midst 
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of the most intense trapped proton radiation belt. As discussed earlier 
in section III.C.3, a power supply other than solar cells (e.g., an RTG) 
would be required. In addition other electronic components would need 
special radiation protection thus incurring some weight penalty. The 
feasibility of this approach would need more detailed study than can be 
given here. 

An alternate approach would be to select a higher equatorial orbit 
outside the most intense radiation region where solar cells would be 
adequately protected. The problem is that as the altitude is increased 
the required power also increases rapidly (figure 29 ) while the rate of 
coverage decreases (figure 18). A reasonable compromise appears to be 
an altitude of about 14,000 km. At this altitude the required average 
power would be about 2.5 kw (figure 29 ) and the satellite would "lap" the 
earth twice daily. A 10 meter antenna could achieve a swath width of about 
280 km at the latitude of the Great Lakes while a 15 meter antenna would 
afford a swath width of about 420 km (figure 28). In either case, however, 
daily coverage of the Great Lakes would not be possible as was the case 
at 6400 km altitude without going to impract ical ly large antennas (greater 
than 30 meters in the range dimension). 

2. Polar (Sun-Synchronous) Orbits 

For polar orbits little is gained by going to high altitudes. The 
coverage rate is low, power required is excessive, and image distortion 
(refs. M , 18) results because the doppler shift produced by the rotation 
of the earth is comparable to the shift produced by the orbital velocity 
of the spacecraft. 
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At low altitudes (400-1000 km) the swath width for single frequency 
operation is on the order of 50-100 km. Rapid coverage of the Great Lakes 
from polar orbit therefore requires mul t i frequency /mul t i beam techniques. 
Typical power requirements as a function of Incidence angle for single 
beam/single frequency operation are shown in figure 32 . 

The curve must be interpreted carefully. The antenna range dimension 
is constantly changing as indicated to suppress the range ambiguity. This 
change is in a direction to increase the antenna gain at low incidence 
angles. Hence, a decline in power requirement results at low incidence 
angles. Also, at low incidence angles larger width radar pulses are 
sufficient to obtain the required range resolution. The attendant band- 
width reduction allows for a further reduction in power. However, these 
effects only delay the inevitable. At 10° incidence the required antenna 
range dimension is 8.5 meters and increases rapidly as the incidence angle 
decreases further. The maximum practical antenna size is quickly encountered 
and no further power reduction is possible. Also the scattering coefficient 
decreases very rapidly at these low angles and eventually i { must dominate 
and cause the power to rise again. 

For polar orbits the swath width falls along the east-west span of 
the Great Lakes. This east-west span is approx imatel y 1300 km. It is 
reasonable to assume that at least four and possibly six beams could be 
synthesized with the type antenna needed for this application. At 20° 
incidence this would require four to six X-band multiple feeds to illuminate 
a 10 x 3 meter antenna with the displacement of the beams along the minor 
axis. 
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Although figure 32 showed the required power to rise rapidly as the 
incidence angle increased, in the range of 10° to 40° this rate is 
relatively slow. In the vicinity of 20°, say ±2° the required power could 
be approximated with an average value. Then superimposing several beams 
over this area and adding the required powers based on the above average 

value will not be much in error. Using this approach, the swath at 20° 

can be increased in multiples of 109 km and 97 watts required power. Four 
beams would then span 436 km with a required 388 watts of power. This 

swath would be sufficient to span the Great Lakes in three days. Six beams 

would imply two day coverage with 582 watts. The approximate required 
power versus coverage period is shown in figure 33. 
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SYSTEM SELECTION 

From the foregoing discussion it appears that daily coverage could 
perhaps be achieved from an equatorial orbit at about 6400 km altitude. 

The 6400 km equatorial orbit "laps" the earth five times daily thus allow- 
ing a daily image to be constructed from a mosaic of five consecutive 
passes. Radar beam swath widths on the order of 250 km are achievable 
at the latitude of the Great Lakes with such an equatorial satellite so 
that an area the size of the Great Lakes (about 900 km north to south) 
could be mapped in five or even four passes. The average power required 
would be about 400 watts. However, because of the intense radiation at 
this altitude (due mostly to trapped protons) a power source other than 
solar cells (e.g., an RTG) would be required. Other electronic components 
would need radiation protection, thus incurring a weight penalty. Such a 
system is judged to be relatively complex and costly and does not appear 
to be an attractive solution. However, this approach would need more 
detailed study than can be given here to adequately assess its feasibility. 

An alternate approach is to select a higher equatorial orbit (above 
11,000 km) where solar cells can be adequately shielded but at the expense 
of daily imaging. An altitude of about 14,000 km appears to be a reasonable 
choice. At this altitude, solar cells can be adequately protected with 
a fused silica glass cover of reasonable thickness. A thickness of about 
20 mils would limit power degradation to 25 percent at the end of five years. 
At 14,000 km the required average power for imaging at the latitude of the 
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Great Lakes would be about 2.3 kw. The satellite would "lap" the earth 
twice daily. At the latitude of the Great Lakes, swath widths of 280-420 
km could be achieved with antenna range dimensions of 10-15 meters 
respectively. Daily coverage of the Great Lakes would not be possible 
since it would require an impractical ly larg4 antenna (greater than 30 
meters in the range dimension). 

Some pertinent parameters for the two equatorial orbits discussed 
above are listed in the following table: 


Equatorial Orbits for Great Lakes Imaging 


Altitude, km 

6410.82 

13916.91 

Period, hours 

3.998 

7.993 

1 ncl i nat ion, deg 

0 

0 

Average power, watts 

400 

2500 

Area scan repetition period , hours 

4.8 

12 

Latitude coverage, deg 

41 -49 

41-49 

Respective beam incidence angles, deg 

26-12 

34-24 

Achievable swath width at 50 deg 
latitude, km 

250 

420 

Antenna size, meters 

12 x 15 

6x15 

Resolution, meters 

100 x 100 

100 x 100 

Receiver effective noise, °K 

500 

500 

Primary power source 

RTG 

Sol ar 
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Conventional s i ngl e-beam/s i ngl e-frequency technology is adequate for 
covering the Great Lakes from an equatorial orbit. However, development 
of a deployable 12 meter x 15 meter or 6 metef x 15 meter antenna would be 
required. This appears to be technically feasible in light of the 
successful development of a IOmeter (approximately) deployable antenna 
for ATS-F. 

To create a mosaic, the radar beam would have to be moved on each 
pass. This could be done by either changing the satellite attitude or 
steering the beam. Steering the beam would be the preferred approach to 
minimize the satellite attitude control requirements. For the small 
steering angle involved (about 6 deg at 6400 kni and 1.5 deg at 14,000 km) 
the beam could probably be electronically steered. Otherwise a slewed 
antenna would be required which would add to the satellite's attitude 
control requirements. 

Eclipsing of the satellite by the earth would occur six times daily 
for the 6400 km orbit and three times daily for the 14,000 km orbit. The 
eclipsing would have no direct effect on the power system design for the 
6400 km orbit since solar cells would not be used. For the 14,000 km orbit, 
eclipsing would occur about nine months of the year with a maximum of about 
10.5 percent of each orbit (about 50 minutes). In terms of earth coverage 
each eclipse would cover about a 25 deg slice in longitude at intervals 
(center to center) of 120 deg. Any one area could be scanned only once 
daily with the satellite in full sunlight, thus requiring three days to 
cover an area the size of the Great Lakes (420 km swath each day). By 



providing battery backup, both of the daily passes over the target area 
could be utilized for imaging, thus halving the time required to 1,5 days. 
However, because of the large amount of eclipsing involved, battery backup 
is impractical. For instance, at least 250 pounds of Ni-Cad batteries 
would be required to provide this backup capability. The batteries would 
also be required to undergo a daily charge/discharge cycle for about nine 
months each year. 

For radar imaging over a range of latitudes from an equatorial orbit 
a slewed antenna would be required. Earth coverage would be restricted 
to latitude bands on either side of the equator. Limiting the radar beam 
incidence angle to a range from 10 to 70 deg, for instance, would limit 
the latitudes to bands from about 10 to 50 deg for the 6400 km orbit and 
about 14 to 62 deg for the 14,000 km orbit. Degradation in performance 
will occur at the lower latitudes because of the change in incidence angle. 
Using the 6400 km orbit as an example, with no change in the radar system 
the swath width would decrease from 250 km at 50 deg latitude to 50 km at 
10 deg latitude. The range resolution would be proportionately degraded 
from 100 meters to 500 meters. Conceivably, this problem could be cir- 
cumvented by using a variable mode radar which would switch to multiple 
beams and shorter pulse widths at the lower latitudes. This would require 
the development of a mu 1 1 i beam/mu 1 1 i frequency synthetic aperture radar which 
may require a substantial technology development effort. As an indication 
of the difficulty involved, the Goodyear Aerospace Corporation has 
estimated that for the EOS (Earth Observatory Satellite) radar, which is 
a single beam/single frequency system, there is only a 50 percent likeli- 
hood of obtaining 5000 hours or about 200 days of operation. The chances 
of obtaining a similar operational capability, much less operations on the 
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the order of several years, with a mu 1 1 i beam/mul t i frequency synthetic 
aperture radar would appear to be quite small at this time. Considerable 
technology development work is indicated in this area. 

An important point to remember is that the required average power is 
a strong function of the target latitude (or beam incidence angle). This 
is illustrated in figure 31 for the 6400 km orbit. Note that while less 
than 400 watts is required at 50 deg latitude about 7 kw would be required 
at 10 deg latitude to utilize the full unambiguous range interval of the 
system. In other words, figure 31 assumes a constant resolution as latitude 
is varied and the antenna range dimension is decreased to match the growth 
in unambiguous range interval. The latitude range covered would thus be 
considerably less than indicated earlier in order to limit the power required 
to a practical value. 

For the sun-synchronous orbit case, the relatively modest power 
requirement allows us to choose an altitude high enough (around 1000 km) 
to completely avoid atmospheric drag and the attendant complexities of a 
drag compensation system. A reasonable choice seems to be the case having 
a three-day coverage pattern and requiring about 400 watts of average power. 

A two-day coverage requiring about 600 watts is also considered practical. 
Daily coverage would require about 1.2 kw which is not as attractive as 
the two-day coverage case in terms of the power required. 

Wide swath coverage from sun-synchronous orbit will also require the 
development of a mul t i beam/mu 1 1 i frequency synthetic aperture radar. The 
number of beams required for the three-day, two-day and daily coverage 
cases would be 4, 6 and 12 respectively. From a practical point of view, 
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12 beams is not consistent with the required antenna range dimension 
of about 3 meters, thus ruling out daily coverage. 

The three-day and two-day coverage cases are summarized below for 
a resolution of 100 meters x 100 meters and a receiver effective noise 
temperature of 500°K. 

Sun-Synchronous Orbits 



Coverage Period 

3-day 

2-day 

Altitude, km 


1011 .5441 

1072.2403 

Period, hours 


1 .75610 

1.77778 

Incl i nat ion, deg 


99.5310 

99.8105 

Average power, watts 


400 

600 

Orbits per 24 hours 


13-2/3 

13-1/2 

Earth track separation 

in longitude, deg 

8.7805 

13.3333 

Earth track separation 

at equator, km 

952.0 

1443.5 

Earth track separation 

km 

at 45 deg latitude, 

655.5 

995.2 

Beams required 


4 

6 

Total achievable swath 

width, km 

436 

654 

Antenna size, length x 

width, meters 

10x3 

10x3 


For sun-synchronous orbits, beam steering is not required since all 
targets eventually come into view resulting in global coverage. This is 
not possible with the equatorial orbits. With global coverage possible, 
more applications for this system would likely result. 
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It should also be noted that the average power levels quoted for 
the sun-synchronous orbits are for a mean incidence angle of 20 deg. 

Unlike the equatorial orbits, the incidence angle can remain nominally 
fixed while still achieving global coverage, thus avoiding the increase 
of power at lower incidence angles that the equatorial orbits experience. 

For a sun-synchronous orbit normal to the sunline at about 1000 km 
altitude, solar eclipsing will occur for a period of about two months 
(roughly mid-May to mid-July) reaching a maximum of about 11 percent of 
the orbit period (about 12 minutes). Battery backup power could be pro- 
vided if imaging is required during this period. About 10-15 pounds of 
Ni-Cad batteries would suffice. Otherwise, batteries would only be 
required for housekeeping functions. 
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CONCLUSIONS 

Applications involving dynamic phenomena such as ice imaging, 
ocean state surveillance, and pollution monitoring require short response 
time (image turnaround time) which in turn implies wide swath width 
coverage. 

Equatorial orbits are not a practical solution. Daily coverage from 
equatorial orbit requires an altitude of about 6400 km which is in the 
heart of the trapped proton radiation belt thus precluding the use of solar 
array power. Orbits above the most intense part of the radiation belt 
such as the case investigated at 14,000 km do not allow daily coverage 
and result in considerably higher power requirements. Battery backup power 
does not appear practical. 

Equatorial orbits do not allow global coverage thus limiting other 
potential appl ications. Coverage of a range of latitudes would require a 
multibeam/multi frequency synthetic aperture radar and a slewed antenna. 

In addition, coverage at the lower latitudes would result in high radar 
beam incidence angles and increased power requirements up to the mul ti- 
ki 1 owatt range . 

Sun-synchronous orbits around 1000 km altitude appear to be the best 
approach to obtaining wide swath coverage. However, a mul t i beam/mul t i frequency 
synthetic aperture radar system will need to be developed. Daily coverage 
of an area the size of the Great Lakes does not appear practical since the 
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large number of beams required (approximately 12) is not consistent with 
the required antenna range dimension of about three meters. Coverage at 
two-day or three-day intervals appears feasible. This would require six 
beams and four beams respectively. The average power required would be 
about 600 and 400 watts respectively. 

Sun-synchronous orbits can obtain global coverage, thus allowing 
potentially more applications than an equatorial orbit. Antenna slewing 
would not be required since all targets eventually come into view. 

Future studies are needed to define this system in more depth. Par- 
ticular attention should be given to assessing the feasibility of a multi- 
beam/mult i frequency synthetic aperture radar and the possibility of using 
a large number of beams such as 12 to allow daily coverage of an area the 
size of the Great Lakes. In addition, studies are needed to identify 
other potential users and to determine what applications are practical. 

The costs and benefits associated with these applications need to be de- 
fined along with the capabilities of alternate approaches such as aircraft, 
in situ observations, etc. to adequately assess the merits of a radar 
imag i ng satel 1 i te. 
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APPENDIX A 

Derivation of V 7 as a Function of Time for a Decaying Orbit 


The angle 4" is defined as the longitudinal relation between the 
sunline and the orbit perpend i cul ar, ^ . Thus, '/ x as a function 

of time can be expressed as 

y = J (o( - ) dt + K (constant) (A-l) 

Since the orbit is taken to be initially sun-synchronous, cL - i 1 for 

« % 

h = h and I - T in equation (1) of the main text. Note that - X^ 
o o a 

As the satellite decays, XI varies with h according to equation (1) 

with I * 1 . 


Thus, 

y = K - J R £ 2 cos I 

Assuming h < < R 


(R + h )~ 3,5 - (R + h )" 3 * 5 


e o 


dt (A-2) 


V = K + 3.5 {p J R e " 2 " 5 cos I o J ( h Q -h) dt 


(A-3) 


The integral of h dt must now be determined. Using h from equation (8) 
of the ma i n text. 


/ 


h dt = 86400 / In 


- 86400 B / 


t \ + <=* & ^o 


v w/c d a , 


dt (A-4) 


where dt is reckoned in days. 
Let 


Z = 


- 86400 B 


(w/c d a ) 


6 h 

+ p* r o 


so that dZ * -86400 B dt 




Substituting into equation (A-4), 
h dt = ( W/C D 


JW/C n A) ( . , n z d2 „ . (W/C d A) 


T* 


~pr 


[z(lnZ-l)j 


(A-5) 
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By substituting equation (A-5) and Z into equation (A-3) and rearranging, 


t'- K = 3.5 iJJl J R e " 2 ' 5 (W/C d A) cos 


In 


-86400B j t \ + 

L lys v 


') cos r I J 

: Li: 

( h j 

e ° 


-8 6400 j t \ 

/* Iw/c d aJ 

“ J S' + 86400 h 


+ ^ o 


° ( M/ W 


(A-6) 


It should be remembered that t is in days and W is in Kgf in equation (A-6) . 
The constant K may be determined by evaluating equation (A-6) at t - 0 
when h = h and ^ . 


\[F J R e " 2 ' 5 (W/C d A) cos I o 

1 

1 

o 

_c 

o 

_c 

> 


L J 


(A— 7) 


L- J 

Combining equations (A-6) and (A-7) we may plot ~ Y versus t 


w/c d a 


W/C d A 


for various initial sun-synchronous orbits (see figure 8). 
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APPENDIX B 


Separation Distance Between Earth Tracks of Two 
Circular Orbits Separted in Longitude by 

Consider a circular orbit about the earth with a given retrograde 
inclination. As the satellite moves along its orbital path the sub- 
satellite point (along radius vector to center of earth) will describe 
a track on the earth's surface. As the earth rotates beneath the 
satellite orbit the resulting successive tracks will move westward across 
the earth's surface. Although the orbit itself describes a great circle 
about the center of the earth the track on the earth's surface will not 
because of the earth's rotation. This fact complicates but does not preclude 
obtaining a closed form approximate solution for the separation distance 
between the tracks. The approach here is to first derive an approximate 
solution based on assuming the tracks are great circles and then refining 
this solution for the actual geometry. 

If we assume that the tracks are great circles, the computation of 
the separation distance involves utilizing a fairly simple set of 
trigonometric formulae. Consider sketch B1 showing two orbit tracks as 
great circles separated in longitude by A/t and having an inclination, T. 
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Let us define the separation distance, S, as measured along a 
great circle on the surface of the earth from the subsatellite point 
at some latitude, X, normal to the current orbit track (orbit track 2), 
to the preceding orbit track (orbit track 1). 

It should be remembered that the geometry in sketch B1 does not 
represent the actual case. In reality the orbit track is continuous as 
the satellite orbits the earth and not a set of discrete circles as is 
depicted. To an observer on the surface of the earth the track would 
appear to cross the equator with an inclination greater than the orbit 
inclination due to the rotational velocity of the earth. However, the 
northernmost position of the orbit track would still reach a latitude con- 
sistent with the orbit inclination, i.e., 'ft -J. The resulting orbit 
tracks thus will not be great circles due to the earth's rotation. However, 
each succeeding track will be skewed or displaced in the same way so 
that for a first approximation we can neglect the earth's rotation in 
calculating separation distance. 











69 


Consideration of sketch B1 and the relationships of spherical 
trigonometry lead to the following series of expressions: 


sin J ~ s i n X 
sin I 


tan d /2 — tan 




cos 


cos I 3 tf ** 

ill 1 ) 

(T - V 


sin ft - sin (T ; \ ^ 0 or ^ = X - X- 0 

sin y 2 

cos "T = sin I cos ( ) 

i 

sin A = sin 1 sin (<P~ + A-U) 
sin r 

cos 9 = - cos T cos ( 7 Y - ft ) + sin Y' sin ( Tf ~ fi ) 

2 2 

cos (X - X ) 

sin d) = sin 'T sin ( X - X ) 

sin 0 

The separation distance, S, can then be found from 

S = R ♦ 


(B-l) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 


(8—6 ) 
(B— 7) 

(B-8) 


where R is the radius of the earth and may be computed for an oblate earth 
e 


from 

R = ab (B-9) 

e .2 .2 .2, .2v, 1/2 

(b +(a - b ) s i n A.) 

where a (equatorial radius) = 6378.160 km 
b (polar radius) = 6356.775 km 
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One way of considering theeffect of earth's rotation on the orbit 
track and thus on the computation of the separation distance is to consider 
the instantaneous direction of motion at any point along the orbit track 
and to imagine that this resulted from an orbit with the appropriate 
inclination about a stationary earth. Thus the orbit track is assumed to 
be a great circle at any instant caused by an appropriately inclined orbit 
and the solution essentially reduces to the one just derived. The problem 
is to compute an instantaneous pseudo- incl i nat ion for any point along the 
orbit track. The longitudinal displacement of the track caused by earth's 
rotation affects all tracks equally and thus does not affect the separation 
d i stance . 

Consider the orbit track as it crosses the equator. As indicated earlier., 
the earth's rotation will cause the track to have a greater inclination 
than the orbit itself. This may be seen by considering the linear velocity 
vector d i agram in sketch B2. 



Sketch B2 

V$ i* the satellite's inertial orbital velocity. V* is the velocity 
of the satellite relative to the earth caused by the relative motion in 
longitude between the earth and the satellite's orbital plane. This 
relative longitudinal motion results from the rotation of the earth and the 
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precession (or regression) of the orbital plane due to oblateness. If 

the angular rotation of the earth is and the nodal precession of 

* 

the orbit is , then the resulting relative angular motion is 

0 

and V* = R 

s 6 

Thus, from sketch B2, the resultant velocity of the satellite relative to 
the surface of the earth Is 

v,. = (V s 2 + V' 2 - 2 v s V 4 ' cos I) 1/2 

and 

V = cos"' (V r 2 + V s 2 - v' 2 ) 

2 V * \ . 

so that at the equator, 

i = r + .v 

For the computation of S at the equator then, equations (B-I) through 

(B-9) developed earlier can be used with I substituted for X. 

For locations other than the equator, imagine that the velocity diagram 

of sketch B2 is simply translated to the subsatellite point at some latitude, 

(refer to sketch Bl). The geometry is basically the same with the 

plane of the velocity vectors normal to the radius vector. Vy lies in the 

direction of the orbit plane and V s is normal to the meridian plane. 

Referring to sketch BJ, the angle between the two would be (at - £) where 

2 

|3 can be solved from the previous series of express ions<> 
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In this case. 


V = 
r 


V + V - 2 V V cos (4+_?/) 
s s s s v 


1/2 - (V 2 H 
s 


V Sin 4 ) 1 


^ is computed as before and from spherical trigonometry. 


» r** - 1 

I = ff - COS 


[sin (S + ^) cos \ 


Here again, I would be substituted for 1 in equations (B- 


v' 2 + 2 V 

s s 


1) through 


(B-9). 
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Figure 1 - Earth, Sun, Satellite Geometry 
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Figure 3 “ Sun-Orbit Geometry for Sunllne on Same Side of Orbit 

Plane as Orbit Perpendicular (Angular Velocity Vector) 
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Figure 8 - Effect of Altitude Decay on (f - V ) as a Function of Time for Various 
Initial Altitudes, u 


h y Km 
o’ 







































































































Prppellant Weight 





Latitude, Degrees 









:Qrb 




























































Figure 20 - Geometry of Distance from Satellite to Earth 
Target Point 
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Figure 21. - Radar images of Lake Erie ice cover on February 22 and 23, 1973. Time of flight over the lake 
was approximately 70 minutes. The center line of the image is the blind ground track of the aircraft. 
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Figure 22 . Photomosaic of ERTS-1 satellite imagery from band 4 showing Lake Erie ice cover on Feb. 17 (eastern basin) 
and Feb. 18, 1973 (central basin). Note the ice movement that occurred in one day along the southern shore. cs _ 68689 



Typical spaceborne radar geometry showing the spacecraft and its ground 
track with the area being imaged to one side of the ground track. AX 
is the apparent separation of the elements of the synthetic array with 
azimuthal dimension, . 




Figure 24 - Geometry showing relationship between incidence angle, © , pulse width, 
T , and range resolution on ground, At . 



Figure 25 - Overall receive pattern for a typical synthetic array. The dotted 
contour is the real antenna pattern and the solid contour is the 
resultant pattern after synthesis. 



Figure 26 - Orbital geometry for SAR showing unambigious 

range interval, AR, the resultant swath width 
and required beam width, SW and /3 , the 
equivalent increment in longitude or latitude, 
for a particular incidence angle © 4 , Q z 
and altitude h. R E is the earth's radius and 
' s the orbital radius. 
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Figure 27 - SAR swath width as a function of SAR 
altitude with SAR velocity as a 
parameter . 
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Figure 28 - Unambigious range interval versus spacecraft altitude 
for X band SAR in equatorial orbit. Target latitude 
48? Maximum interval was fixed at 900 km to span the 
Great Lakes. Also shown are the number of passes 
daily over this target. 



Required Average Power , 



5000 1 0 , 000 20,000 50, 000 

Spacecraft Altitude, Km 


Figure 29 - Required average power versus orbit altitude for 
X band SAR in equatorial orbit. Target latitude 
48 °. 
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Figure 30 - Required local oscillator stability versus spacecraft 
altitude for X band SAR in equatorial orbit. Target 
1 at i tude 48° . 



Required Average Power, Kw 


10.0 


1 .0 


~ 

- 

— 1. - - l L : ! I J 

0 10 20 30 40 50 60 

Target Latitude, Degrees 

Figure 31 - Required average power to utilize full 
unambigious range interval for X band 
SAR at 6400 Km equatorial orbit versus 
target 1 at i tude . 
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Figure 32 - Required average power and antenna range 

dimension for X band SAR In 1000 km altitude 
polar orbit as a function of incidence angle. 
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Figure 33 ~ Required average power for X band SAR 
in 1000 km polar orbit versus coverage 
period. Powers are approximate only. 
Mean incidence angle 20°. 


